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Abstract. In fractional calculus there are two approaches to obtain fractional derivatives. The first 
approach is by iterating the integral and then defining a fractional order by using Cauchy formula 
to obtain Riemann fractional integrals and derivatives. The second approach is by iterating the 
derivative and then defining a fractional order by making use of the binomial theorem to obtain 
Griinwald-Litnikov fractional derivatives. In this article we formulate the delta and nabla discrete 
versions for left and right fractional integrals and derivatives representing the second approach. 
Then, we use the discrete version of the Q-operator and some discrete fractional dual identities to 
prove that the presented fractional differences and sums coincide with the discrete Riemann ones 
describing the first approach. 
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1 Introduction and preliminaries 

Fractional calculus (FC) is developing very fast in both theoretical and applied aspects. As a result 
FC is used intensively and successfully in the last few decades to describe the anomalous processes 
which appear in complex systems (T] [2] |3] E] El |6] . The complexity of the real world phenomena is a 
great source of inspiration for the researchers to invent new fractional tools which will be able to dig 
much dipper into the mysteries of the mother nature. Historically the FC passed thorough different 
periods of evolutions and it started to face very recently a new provocation: how to formulate properly 
its discrete counterpart [g|8]|9l[l0][IT][T3[l3l[ll[t5][l6][I3[l8l[ig At this stage we have 

to stress on the fact that in the classical discrete equations their roots are based on the functional 
difference equations, therefore the natural question is to find the generalization of these equations 
to the fractional case. In other words we will end up with generalizations of the basic operators 
occurring in standard difference equations. As it was expected there were several attempts to do this 
generalization as well as to apply this new techniques to investigate the dynamics of some complex 
processes. In recent years the discrete counterpart of the fractional Riemann-Liouville, Caputo, were 
investigated mainly thinking how to apply techniques from the time scales calculus to the expressions 
of the fractional operators. Despite of the beauty of the obtained results one simple question arises: 
can we obtain the same results from a new point of view which is more simpler and more intuitive? 
Having all above mentioned thinks in mind we are going to use the binomial theorem in order to get 
Griinwald-Litnikov fractional derivatives. After that we proved that the results obtained coincide 
with the ones obtained by the discretization of the Riemann-Liouville operator. In this manner we 
believe that it becomes more clear what the fractional difference equations bring new in description 
of the related complex phenomena described. 

For a natural number ra, the fractional polynomial is defined by, 

t(") = ff (t - j) = ^^(* + ^^ ^ , (1) 

where F denotes the special gamma function and the product is zero when t -|- 1 — j = for some j . 
More generally, for arbitrary a, define 

^(* + ^^ , (2) 
r(t + 1 - a) ' ^ ' 

where the convention that division at pole yields zero. Given that the forward and backward differ- 
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ence operators are defined by 

A/(i) = fit + 1) - fit) . SJfit) = fit] - fit - 1) (3) 

respectively, we define iteratively the operators A™ = A(A'"~-'^) and V" = V{V"'~^), wliere m is 
a natural number. 

Here are some properties of the factorial function. 

Lemma 1.1. fW}l) Assume the following factorial functions are well defined, 
(i) At(") = 

(a) it - ^i)t'-'^^ = t(f+i), where ^ G M. 

(Hi) = ^i^l + l'). 

(iv) If t <r, then t'"' < r'"' for any a> r. 

(v) IfO<a< 1, then > 

(vi) = (i-/3)(")t(/3). 

Also, for our purposes we list down the following two properties, the proofs of which arc straight- 
forward. 

Vsis - = (a - l)(p(s) - (4) 

Vt(p{.) - = -(a - \)ipis) - t)("-^l (5) 

For the sake of the nabla fractional calculus we have the following definition 

Definition 1.1. (IM \M\M\M) 

(i) For a natural number m, the m rising ( ascending ) factorial of t is defined by 

•m — l 

= ]^ it + k), t^=l. (6) 

fc=0 

(ii) For any real number the a rising function is defined by 

^- ^ r(f + a) {.. 2,-1,0}, 0"=0 (7) 

Fit) 

Regarding the rising factorial function we observe the following: 

(i) _ 

V(t") = (8) 

(ii) 

(t") = (t + a-l)('^). (9) 

(iii) _ 

Atis - pit)r = -ais - pit)r-^ (10) 

Notation: 

(i) For a real o > 0, wc set n = [a] + 1, where [a] is the greatest integer loss than a. 
in) For real numbers a and b, we denote Na = {a, a + 1, ...} and (,N = {^i ^ ^ 1> ■■■}• 
iiii) For n G N and real a, we denote 

eA'V(t) ^ (-l)"A"/(t). 

iiv) For n G N and real b, we denote 

vg/W ^ (-i)"v"/(i). 

Definition 1.2. Let ait) = t + 1 and pit) = t — 1 be the forward and backward jumping operators, 
respectively. Then 

(i) The (delta) left fractional sum of order o > (starting from a) is defined by: 

1 

A-"/{t) = — teN,+,. (11) 



(ii) The (delta) right fractional sum of order a > (ending at b) is defined by: 

b b 



6 6 

i,A-"/(i) = — {^-'^{*))'"""/W = ^ E (pW-i)<"-''/W, (12) 



.i^^. r(a) ^^^^^ 
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(iii) The (nabla) left fractional sum of order a > (starting from a) is defined by: 

v-"/(i) = ^ iZ teN„+i. (13) 

(iv) The (nabla) right fractional sum of order o > (ending at b) is defined by: 
6-1 6-1 

Regarding the delta left fractional sum we observe the following: 

i) Aa"' maps functions defined on Na to functions defined on Na+a- 

ii) u{t) = Aa" f{t), n g N, satisfies the initial value problem 

A"-u{t) = f{t), t e Na, u{a + j ~1) =0, j = l,2,...,n. (15) 

iii) The Cauchy function — — vanishes at s = t — (n — 1), ...,t — 1. 

Regarding the delta right fractional sum we observe the following: 

i) bA~°' maps functions defined on to functions defined on 6-aI^- 

ii) u{t) = f,A~"/(t), n G N, satisfies the initial value problem 

Vg«(t) = /(t), t€ tN, uib-j + l) = 0, j = l,2,...,n. (16) 

iii) the Cauchy function vanishes at s = t + 1, t + 2, t + (n — 1). 

Regarding the nabla left fractional sum we observe the following: 

i) Va " maps functions defined on Nq to functions defined on N^. 

ii) Va "/(t) satisfies the n-th order discrete initial value problem 

V"2/(t) = /(t), V'y(a)=0, i = 0, 1, n - 1 (17) 

iii) The Cauchy function '*~r'(n)" ^ satisfies V"{/(i) = 0. 

Regarding the nabla right fractional sum we observe the following: 

i) 6^~" maps functions defined on j,N to functions defined on j,N. 

ii) b^~"f{t) satisfies the n-th order discrete initial value problem 

eA"y(t) = /(t), eA'y{b) = 0, i = 0,l,...,n-l. (18) 
The proof can be done inductively. Namely, assuming it is true for n, we have 

eA"+i 6V-("+i)/(t) = eA"[-A 6V-("+i'/(t)]. (19) 
By the help of JTOj, it follows that 

eA"+i 6V-("+i)/(t) = eA" 6V-"/(i) = fit). (20) 

The other part is clear by using the convention that 5I]fc=s+i = 0- 
(iii) The Cauchy function ^^~r\n)"~^ satisfies eA"j/(t) = 0. 

Definition 1.3. (i)f8^ The (delta) left fractional difference of order a > (starting from a ) is 
defined by: 



t — (n — (x) 

r(n - a) 

(ii) fl5( The (delta) right fractional difference of order o > (ending at b ) is defined by. 



Ag/(t) = A"A-'"-"'/(t) = _ , E {t-^(s)Y""''^^f(s), iGN„+(„_,) (21) 



5A"/(i) = V§ 6A-("-)/(i) = J2 («-^W)<"-"-"/(^), t6 6-(n-c)N (22) 

(iii) \16\l The (nabla) left fractional difference of order a > (starting from a ) is defined by: 

V^/(t) = V"V-("-°V(t) = ^7 r 5^ (t-p(s))"-°-V(«), teN„+i (23) 

r(n-a) ^i^^ 
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(iv) f \25f . I27j The (nabla) right fractional difference of order a > (ending at b ) is defined 

by: 



b-l 



i,V"/W = qA" 6V-("-"'/(i) = ^-^ J2(^-p(t)r-'-~'f(s), iG ,_iN (24) 

Regarding the domains of the fractional type differences we observe: 

(i) The delta left fractional difference A" maps functions defined on Na to functions defined on 



(ii) The delta right fractional difference i,A'* maps functions defined on [,N to functions defined 



b— (n — a) 



N. 



(iii) The nabla left fractional difference V° maps functions defined on Na to functions defined 

on Na + n ■ 

(iv) The nabla right fractional difference V" maps functions defined on [,N to functions defined 
on b_„N . 

Lemma 1.2. Let 0<7i— l<a<n and let y{t) be defined on ^a■ Then the following 

statements are valid. 

(i) {K)y{t - a) = "^a-iVii) for t e N„+,. 

(ii) {A-")y{t + q) = V-°iy(t) for t g N^. 

Lemma 1.3. \2T^ Let y(t) be defined on t+iN. Then the following statements are valid. 

(i) ( bA")yit + a) = 6+1 V°y(t) for t e i,_„N. 

(ii) ( bA-")y{t -a)= 6+iV-"j/(t) for t 6 bN. 

If f{s) is defined on Na fl i,N and a = b {mod 1) then (Q/)(s) = /(a + b — s). The Q-operator 
generates a dual identity by which the left type and the right type fractional sums and differences 
are related. Using the change of variable u = a + b — s, in |14) it was shown that 

A-"Q/(t) = Q sA— /(t), (25) 

and hence 

AZQf{t) = (Q t,A-/)(i). (26) 
The proof of ((26j follows by the definition, lj25p and by noting that 

-QV/(t) = AQ/(t). 
Similarly, in the nabla case we have 

V-"Q/(t) = Q6V-"/{t), (27) 

and hence 

KQfit) = (Q 6V"/)(i). (28) 
The proof of II28I I follows by the definition, II27I I and that 

-QAfit) = \7Qf(t). 

For more details about the discrete version of the Q-operator we refer to [271 . 

From the difference calculus or time scale calculus, for a natural n and a sequence /, we recall 

n 

A"f(t) = J2i-^)H I )f{t + n-k), (29) 

fc=0 

and 

n 

V"/(t) = 5^(-l)*( ^ )/(t-n + fc), (30) 

fe=0 

2 The fractional differences and sums with binomial 
coefficients 

We first give the definition of fractional order of II29I I and II30I I in the left and right sense. 

Definition 2.1. The (binomial) delta left fractional difference and sum of order « > for a function 
f on defined on Na, are defined by 
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a + t — a 

a-\-t — a 
fc=0 

a+k -1 , 

Definition 2.2. Tfce (binomial) nabla left fractional difference and sum of order o > for a 
function f on defined on N^, are defined by 

• H 

t-a — l 
k-0 

• (b) 



Analogously, in the right case we can define: 

Definition 2.3. The (binomial) delta right fractional difference and sum of order o > for a 
function f defined on jN, are defined by 

• H 

a + b-t 

bAB" = ^ (-!)'=( ^ )/(t-a + fe), t&b-n + aN 
fc=0 

• (b) 

-a + b-t 

bAB-°'= E )/(* + " + fc), teb-aN 

k = 

Definition 2.4. The (binomial) nabla right fractional difference and sum of order o > for a 
function f on defined on 5N, are defined by 

• H 

b-t-l 

bVB'-= E (-!)'( t )fit-k), ie 5_„N 

fc=0 

• (b) 

b-t-l 

,vs-"= E ~k )/(*-^)' *e 

fc=0 

We next proceed to show that the Riemann fractional differences and sums coincide with the 
binomial ones defined above. We will use the dual identities in Lemma 1 1.2 1 and Lemma ll.3l and the 
action of the discrete version of the Q-operator to follow easy proofs and verifications. In 1161 the 
author used a delta Leibniz's Rule to obtain the following alternative definition for Riemann delta 
left fractional differences: 

t + c, 

A^/(i) = - -E(*-<^W)'""""' tef^a + n-a., (31) 

then proceeded with long calculations and showed, actually, that 

A2f(t) = BA^f(t), A--/(t) = iJA-"/{t) (32) 

Theorem 2.1. Let f be defined on suitable domains and a > 0. Then 

• 1) 

A^f(t) = BAZf{t), A-"/{t) = BA--/(t) 



• (b) 

where (— 1)'°( ? ) = 
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• 2) 

6A"/{t)= 5AB"/{t), 5A-"/(t) = AB-«/(t) 

• 3) 

• 4) 

bV°'fit)= tVB°'fit), ^V~°'fit)= t^B-''f{t) 

Proof. • 1) follows by ll32l l. 

• 2)By the discrete Q-operator action we have 

CK + f — a 

6A-/(t) = QA,{Q/)(i) = Q Yl t ){Qmt + a-k)= tAB'^fit). 

k = 

The fractional sum part is also done in a similar way by using the Q-operator. 

• 3) By the dual identity in Lemma 11.21 (i) and II32I I , we have 

^2 fit) = A^+i/{t + a) = BA2+^f{t + a) = BVZf{t)- 
The fractional sum part can be proved similarly by using Lemma ll.2l iii) and I I32I I. 

• 4) The proof can be achieved by either 2) and Lemma 11.31 or alternatively, by 3) and the 
discrete Q-operator. 

□ 

Remark 2.1. In analogous to \31\l the authors in I2ffj used a nabla Leibniz's Rule to prove that 

yaf{t) = ^^, it-pi')f^fi')- (33) 

In \25f the authors used a delta Leibniz's Rule to prove the following formula for nabla right 
fractional differences 

1 6-1 

Similarly, we can use a nabla Leibniz's Rule to prove the following formula for the delta right 
fractional differences: 

1 * 

We here remark that the proofs of the last three parts of Theorem 12. 1\ can be done alternatively 
by proceeding as in \16fl starting from i33\ }, (34^ and iSSV , Also it is worth mentioning that mixing 
both delta and nabla operators in defining delta and nabla right Riemann fractional differences was 
essential in proceeding, through the dual identities and the discrete Q-operator or delta and nabla 
type Leibniz's Rules, to obtain the main results in this article JgTf . 
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